Take Off Environment OEI AFM Procedure

Introduction

This notebook deals with a realistic model of the aircraft at take-off conditions. We deal with transport-type airplanes and write the general take-off equations. There are several different types of aircraft take-off. We will refer exclusively to Conventional Take-Off.

Considerations

We consider a take-off operation on a vertical plane. Due to the level of approximations that it is possible to achieve, in the calculation of the take-off of a conventional airplane, we will assume that the runway reference system and the ground reference system are equivalent.

The sequence of airplane speeds during take-off are given in the figure below.

img/takeoff_phases.png

The airplane accelerates, the nose landing gear lifts off and the airplane rotates with a speed $V_R$; the airplane lifts off the ground with a speed $V_{LOF}$ and follows a rectilinear flight path to clear a conventional screen at 35 feet from the ground.

The take-off is divided into three segments:

Ground Run

The ground run is the distance between brake-release and lift-off of the forward wheels. To calculate the ground run, we write the dynamics equations on the centre of gravity of the aircraft in the horizontal and vertical directions, for a take-off from a horizontal runway. The engine thrust is aligned with the vector velocity.

Governing equations:

State Increments:

The thrust is given by the thrust model which takes as an input the speed.

Target: Reach $V_R$

Rotation and Initial Climb

The lift-off point is reached with at least one landing gear on the ground. The rotation of the aircraft must be done with a small angle or else there is a risk of a tail strike on the runway. The rotation run is the distance between lift-off of the forward wheels and lift-off of the main landing gears.

Governing equations:

Constant: $q = Rotation Rate$

State Increments:

Target : $W = L(\alpha)$

Airborne phase to the height of the screen ($V_{LOF}$ up to $V_{35ft}$)

img/climb.png

In this segment we will try to compute the state of the aircraft from $V_{LOF}$ up to point A.

Governing equations:

Target: Reach 35ft

Remark: In this model the incidence angle $\epsilon$ is considered to be 0.

The climb velocity in general accelerated flight is

$v_c = \frac{T-D}{W}U - \frac{U}{g}dU/dt$

Our goal is to ascend with a max rate of climb, which means minimizing the acceleration. The pilot can manage the acceleration of the aircraft by adjusting the stick input, which alters the change of $\alpha$. A simple gain control loop is the most straightforward method to simulate this process.

A simple schema of the control function is illustrated below:

img/control_loop.png

TakeOff Procedure

Import Libs

Inputs

Methodology Part 1

Before initiating the TakeOff Procedure some preliminary steps have to be made first.

A dedicated class was created to prepare the take-off

Calculate Stall Speed

We calculate the stall speed with the following equation: $ V_{stall} = \sqrt(\frac{2 * mg}{(\rho S C_{Lmax}}) $

Define the rest of the characteristic speeds

Check that we are compliant with the JAR25.121b Regulation

In the case of non compliance of $V_{2min}$ with JAR25.121b, we take the ssg of JAR(2.4%) and we calculate the new $V_2$

All these methods are summed up into one function in the same class called pilot_preparation. pilot_prep-3.png

Methodology Part 2

With our characteristic speeds established, we can commence the takeoff process. The initial step is to start with the minimum $V_R$ and gradually attain $V_2$ while ascending to a height of 35ft.

If any of the flight's corresponding speeds fail to meet the criteria, we increase $V_R$ by $1\% V_{stall}$ and initiate the iterative process anew until all the criteria are satisfied.

This class inherits the three takeoff segments, each of which inherits the takeoff preparation class. If the $V_R$ fails to meet the takeoff criteria, we increment it by 1% of the stall speed until all the criteria are met.

Create the object and run the simulation

PostProcessing

So, the algorithm finally converged to a specific $V_R$!

In order to post process the results we are going to transform the log dictionary to a pandas dataframe.

Characteristic Instants

Now we are going to check the time and speed values that the algorithm found about the characteristic speeds.

We can see that the rotation started 30 seconds after the break release.

Then 4 seconds later the aircraft lifts off!

It only took about 9 seconds for the aircraft to reach an altitude of 35ft!

We notice that $V_{LOF}$ is very close to $V_2$, which makes sense since the goal on the airborne phase is to utilize most of the aircraft's kinetic energy for climbing, resulting in a steep 𝑉𝑧(𝑡) curve.

Remark: There is a distinction between the values $(V_2, V_{35ft}) $. This discrepancy is caused by the controller's impact. The acceleration never reaches zero, but instead oscillates around that value. Consequently, we can anticipate slight differences in the speed values during the airborne phase.

Initial characteristic speeds

Remark: The $V_R$ is not exactly the same with the above results due to the timestep's sensitivity. If we reduce the timestep then the increments dx and dv should tend towards continuity, resulting in a more precise value.

Plots

$x = f(t)$

The above curve has the anticipated shape, with the aircraft following a non-constant accelerating motion during the ground segment and a mostly linear motion during the airborne phase, where zero acceleration is the objective.

Specifically, the curve displays a concave-up shape up to 34 seconds, after which it tends to become more linear.

$TAS = f(t)$

The above graph depicts the velocity increasing in a non-linear fashion up to 34.2 seconds, which is expected since the acceleration is not constant during the ground phase. Subsequently, it stabilizes at the 𝑉2 value, which is logical given that the pilot aims to attain zero acceleration during the airborne phase to facilitate the fastest possible ascent.

$V_z = f(t)$

As expected the rate of climb is 0 until lift-off(34.2 secs) and then we see a rapid increase due to the maximazation of SEP(Specific Excess Power). The peak of this curve marks the peak of the climb gradient, since $V_z = TAS\sin\gamma$ and TAS is quasi static.

Angles Graph

In this graph we can see that until lift-off $\alpha$ and $\theta$ have the same value since $\gamma$ is 0. When the a/c is airborne and tries to climb the value of gamma goes up. $\gamma$ is approximately given by the following formula:

$\gamma = \frac{T\cos\alpha -D}{W}$

This formula makes evident that in order to increase $\gamma$ we have to lower $\alpha$. This happens due to two main reasons.

For $\alpha$ we see that it increases until it reaches the target acceleration and then starts to decrease.

$\theta$ is calculated by the known formula :

$\theta = \alpha + \gamma$

$ \alpha, aoa = f(t)$

In the airborne segment, the pilot uses 𝛼 to control the acceleration of the aircraft. When the aircraft is airborne with a positive acceleration, the pilot pulls the stick to increase 𝛼 and decrease the acceleration. The peak of 𝛼 is reached when the acceleration is close to zero. Subsequently, the acceleration becomes negative, and the pilot pushes the stick to decrease 𝛼 and control the acceleration close to zero.

Conclusion

This notebook presents the AFM OEI Take Off procedure. The user has to input the initial a/c and atmospheric conditions and then the code is able to produce an analytical decription of the take off procedure. The results are satisfying and produce an accurate description of the take off procedure. Then the user can post process the information as he likes.

Limitations

In this simulation there are several factors that were not taken into account:

These factors limit the accuracy of the model

References